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ABSTRACT 

It is believed that the biasing of the galaxies plays an important role for 
understanding the large-scale structure of the universe. In general, the biasing of 
galaxy formation could be stochastic. Furthermore, the future galaxy survey might 
allow us to explore the time evolution of the galaxy distribution. In this paper, the 
analytic study of the galaxy-mass density relation and its time evolution is presented 
within the framework of the stochastic biasing. In the weakly non-linear regime, 
we derive a general formula for the galaxy-mass density relation as a conditional 
mean using the Edgeworth expansion. The resulting expression contains the joint 
moments of the total mass and galaxy distributions. Using the perturbation theory, 
we investigate the time evolution of the joint moments and examine the influence of 
the initial stochasticity on the galaxy-mass density relation. The analysis shows that 
the galaxy-mass density relation could be well-approximated by the linear relation. 
Compared with the skewness of the galaxy distribution, we find that the estimation 
of the higher order moments using the conditional mean could be affected by the 
stochasticity. Therefore, the galaxy-mass density relation as a conditional mean should 
be used with a caution as a tool for estimating the skewness and the kurtosis. 



Subject headings: cosmology:theory— large scale structure of universe, biasing 



1. Introduction 

The enormous data of the large scale galaxy distribution will be obtained in the future. It is 
expected that the physics in the early universe can be explored by using these data. According to 
the standard picture based on the gravitational instability, the primordial fluctuations produced 
during the inflationary stage have evolved into the large scale structure observed in the galaxy sky 
survey. To compare the theoretical prediction for the density fluctuation with the observation, we 
need the relation between the total mass and the galaxies. Due to the lack of our knowledge of the 
galaxy formation process, however, the statistical uncertainty between the galaxies and the total 
mass density arises. This uncertainty affects the determination of the cosmological parameters 
from the observation ( |Hamilton 1997 ). 
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Denoting the fluctuation of total mass as 5 m and that of galaxy distribution as 8g 1 the linear 
biasing relation 

$g = b5 m , (1) 

is frequently used in the literature. Here, b is the biasing parameter. A naive extension of ([!]) to 
the quasi non- linear regime is obtained by taking the expansion in powers of 8 m : 

5 9 = f(S m ) = J2 b ^(6m) n . (2) 

n=l 



The non-linear biasing parameters b n is observationally determined by several authors ( [Fry 



Gaztahaga 1993 , Fry 1994 , |Gaztahaga Frieman 1994 ). However, we must note that there is an 



assumption behind equations ((T|) and @ that the statistics of the galaxy distribution is completely 
characterized by that of the total mass S m . In reality, the galaxy formation process could be 



stochastic ( Cen fc Qstriker 1992 ). Stochastic biasing is a notion to treat this situation. Note 
that the deterministic biasing corresponds to a special case of the stochastic biasing. Stochastic 
property of the galaxy biasing is recently studied by many authors. Pen (1997) has discussed an 
observational method to determine the stochasticity. Scherrer &; Weinberg (1998) has studied the 
influences of the non-linearity and the stochasticity on the statistics for the two point correlation 
function. 

Due to the gravitational dynamics and galaxy formation process, the relation of the galaxies 
and the total mass is generally non-linear and time dependent. The recent observation at high 
redshift shows that the galaxies at z — 3 are strongly biased. The evidence of the strong biasing 



at high redshift suggests that the time evolution of the biasing should be considered ( Peacock 



et al. 1998 ). Furthermore, Magliocchetti et al. (1998) showed that the measurement of the 
angular correlation function of the radio galaxies obtained from the Faint Images of the Radio 
Sky at Twenty centimeters is consistently explained by taking into account the evolution of the 
biasing. The time evolution of the galaxy biasing has been firstly studied by Fry (1996). Tegmark 
& Peebles (1998) extended it to the stochastic biasing and examined the linear evolution. In 



the previous paper ( [Taruya, Koyama fc Soda 1999 ), we have investigated the quasi non-linear 



evolution of the stochastic biasing. 

A general formalism for the stochastic biasing is recently proposed by Dekel & Lahav (1998). 
The non-linearity and the stochasticity of the galaxy biasing are measured by the conditional 
mean and its scatter. They performed the numerical simulation and evaluated the conditional 
mean of 5 g at a given total mass fluctuation 5 m . Blanton et al. (1998) studied the galaxy-mass 
density relation using the hydro dynamical simulation and explore the physical origin of the scale 
dependent biasing on small scales. 

The purpose of this paper is to investigate the galaxy-mass density relation and its time 
evolution analytically. Using the Edgeworth expansion, we will derive a general formula for the 
conditional mean. Our main result is equation (|25|). The final expression can be written in terms 
of the joint moment of the galaxy and the total mass. Since the gravity induces the time evolution 
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of the joint moments, the galaxy-mass density relation cannot be static. Neglecting the merging 
and galaxy formation process, we calculate the joint moments using the perturbation theory of 
the Newtonian cosmology. Then we study the time evolution of the galaxy-mass density relation. 
We find that the estimation of higher order moments such as skewness using the conditional mean 
could be affected by the stochasticity. 

We organize this paper as follows. In section ||, the basic description of the stochastic biasing 
is explained. Within the formalism, a general formula for the galaxy-mass density relation is 
derived. Then we examine the time evolution of the galaxy-mass density relation in section ||. The 
final section is devoted to the summary and the discussions. 



2. Stochastic biasing and galaxy-mass density relation 
2.1. Basic description 

Let us recall the definitions of 6 m and 6 g : 

S m (x) = S 9 ( X ) . VMz^L, (3) 

p n g 

where variables with overbars indicate the homogeneous averaged density. In this paper, we treat 
the smoothed density fields 5 m (R) and S g (R) using the spherical top-hat window function Wr(x) 
with the smoothing radius R: 



S m>g (R) = J d 3 x W R (x)5 my9 (x). (4) 

Importantly, the one-point distribution function itself depends on the smoothing scale. As the 
smoothing scale becomes larger, the distribution approaches to the Gaussian distribution. 

The statistical information of the galaxy biasing is obtained from the joint probability 
distribution function (PDF) of the galaxy and the total mass distribution. In our formulation, 
5 m and 5 g are regarded as the independent stochastic variables. To represent their stochasticity, 
we introduce the auxiliary random fields A m (x) and A g (x). On large scales, the smoothed fields 
8 m ,g(R) could be regarded as nearly Gaussian fields. The density fields can be expressed as 
functions of the Gaussian variables and expanded in powers of these variables whose variances are 
small enough. Denoting such Gaussian variables as A m ^ g , the distributions 5 m and 5 g are given by 

S m = f(A m ,A g ), 5 g =g(A m ,A g ). (5) 

Here, we wrote equation (||) in somewhat schematical way. The "functions" / and g are not 
assumed as the local functions. The generating function for the one-point functions S m (R) and 
5 g (R) becomes ( [Taruya, Koyama fc Soda 1999 ) 

Z(J m ,J g ) = ^ e i{JmSm{R)+J a S g {R)}\ 
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J J VA m {x)VA g (x) AA _1 exp 



-(A m ,A 9 )G" 1 



A„ 



x exp 



d 3 XW R (x) {J m 5 m (A m! g) + Jg5g(A mj g)} 



(6) 



where N is a normalization constant and G is a 2 x 2 matrix which gives the stochastic property 
of A m and A g . The joint moments of 5 m and 5 9 are deduced from the generating function as 



([5 m (RW[S 9 (R)} k ) = r^ 



dJ° m dJ k g 



(7) 



The ensemble average (• • •) is taken with respect to A m and A g . Once we obtain the joint moments 
(|7|), the joint PDF V(5 m , 5 g ) for the smoothed density fields 6 m (R) and 5 g (R) is constructed: 



V(8 m ,Sg 



(2tt)= 



dJ rri dJg Z^Jm-iJg) & 



(8) 



where we simply denotes 8 m ,g{R) as S mt 



a- 



We should keep in mind that the gravitational instability usually affects the galaxy and the 
total mass distribution. Accordingly, the functions / and g in the assumption (g) are determined 
by the gravitational dynamics and they would be generally non-local and time dependent (see 
section ||). Since the resulting joint PDF also depends on time, we would have the time evolving 
biasing relation. 



2.2. Stochastic biasing in the weakly non-linear regime: conditional mean and 

biasing scatter 

To proceed further, we shall employ the perturbative approach, hereafter. Taking the 
variables A m)S as seeds of perturbation, S g and S m can be expressed as 

5 m! g(A m ,Ag) = 5$g + 5%g + ---, (9) 

where the n-th order perturbed quantities Smjg are of the same order as [A m>g ] n . Substituting the 
expansion @ into the joint moment (0), the non- vanishing lowest order quantities become the 
second moments. They are characterized by the three parameters: 

/rjf(i)i2\ ttW 1 ^ 

(t^l 2 ) v / ([4 1) ] 2 )([^ ) ] 2 ) 

cr m is regarded as the variance of the total mass, b is the linear biasing parameter and r is the 
cross correlation. Note that the parameter r reflects the stochastic property of 5 m and 6 g . At the 
lowest order level, the deterministic biasing relation (|l]) holds if r = 1. 
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The higher order perturbations leads to the non-linear galaxy biasing which describes the 
relation between the non-Gaussian distributions of galaxy and the total mass. Note that the 
higher order correction also affects the stochastic property of 5 m and 5 g . To investigate the galaxy 
biasing in the weakly non-linear regime, we should analyze the non-linearity and the stochasticity 
separately. 

Several authors proposed to use the conditional mean of the galaxy distribution defined by 
flDekel fc Lahav 1998| , |Blanton et al. 1998|) 



(6 9 }\5 m = / dS g 6 g V(6 g \6 m ) ; V{5 g \8 m ) = rK Jl ™ J , (11) 



T ? {8g, S m j 

J g vg>'\vg\vm) , ryu g \v m) - — 

where V(5 g \5 m ) denotes the conditional PDF at a given 5 m . From the conditional mean {6 g )\s m , 
we can know the non- linear relation of the biasing as the function of 5 m . We should keep in mind 
that the conditional PDF has a width around the conditional mean ($ g )\s m - The biasing scatter 
can be introduced to measure the stochasticity in the biasing relation ( pekel &: Lahav 1998 ): 



^ = S g -(S g )\ Sm . (12) 

The stochastic property of the non-linear biasing relation is understood from the variance and the 
higher moments of e. 

Consider the density fields on large scales. Owing to the assumption (||) and the perturbative 
treatment (||), we can obtain the analytic expression for the conditional mean (5 g )s m - The 
non-linear galaxy-mass density relation using the Edgeworth expansion leads to the expression 

(S g )\s m = c 1 S m + ^(S 2 m -(S 2 m }) + ---. (13) 

We will show in the next subsection that the coefficients c\ and C2 are expressed in terms of the 
joint moments of 8 m and 5 g (see eq.p6|). 

Before deriving the galaxy-mass density relation, we mention the differences between the 
relation ( |l3| ) and (g). From the relation (||), Fry Sz Gaztahaga (1993) derived the useful formulae 
which relate the higher order moments of the galaxies with that of the total mass. For the third 
moments, we have 

S3, 9 = l(S 3 , m + 3^), (14) 

where the linear biasing parameter b\ is simply denoted as b. The quantities Ss ; m and S^g are the 
skewness of the total mass and the galaxy distribution, respectively. In the tree level perturbation 
theory, these can be written as 



For the galaxy-mass density relation (fl3D, we can also obtain the equation similar to ( |14D (see 
eqs.[^6| and [p8| ). We will see in section that the formula ( |l4|) is recovered when the cross 
correlation r becomes unity. Hence, the galaxy-mass density relation (|l~3|) can be regarded as a 
generalization of the relation (|2|) examined by Fry & Gaztahaga (1993). 
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2.3. Derivation of galaxy-mass density relation: a general formula 

We are in a position to derive the galaxy-mass density relation (Jl3|) as a conditional mean. 



In section 2A, we defined the generating function Z(J m , J g ) to construct the joint PDF V(S m , 8 g ). 
Using Z, we obtain the generating function for the connected diagrams: 

W(J m ,J g )^lnZ = £ (iJrn^fmy (iJ 9 ^)) k , 



where we define 



jk {{5m) 2 y 2 <(<y 2 ) fc / 2 ' 

(• • -) c denotes the connected part of the moments. By definition, we have A02 = A20 = 1- The 
perturbative treatment @ implies that Xij satisfies the following scaling relation. Using the 
variable a m defined by equation ( [l0[ ) and assuming that the lowest order variance of the galaxies 
([4 1} ] 2 ) is of the same order of magnitude as that of the total mass, we have 

\ jk = S jk o£ k - 2 + 0(o£ k ). (16) 

The variable Sij is of the order of unity and we obtain S20 = <So2 = 1 ■ 

As long as the variance a m is small, treating a m as an expansion parameter, we can expand 
the generating function W. Thus we proceed to approximate the joint PDF V(5 g ,5 m ). The 
expansion of W is referred to as the Edgeworth expansion ( Juszkiewicz et a/.1995| , [Bernardeau & 
Kofman 1995| , Chodorowski &: Lokas 1997| ). Let us write the conditional mean (6 g )\s m as 



<Vk = ; = Jd5 g 5 g V(8 m , S g ). (17) 



V(5, 

Now, define the normalized variables: 



V = J m \[^L), u = J g ^), x=-^=, y =-?£=. (18) 

V™ \ \ 5 q) 



Using the generating function W, the numerator and the denominator of equation ( |l7| ) are 
expressed in terms of the above variables: 



V{S m ) = — ( t^tV 72 / dfie-^ x+w ^'°\ (20) 



A 



27T \(Sl 

Applying the Edgeworth expansion to the above equations, we can obtain the galaxy-mass density 
relation. The details of the calculation are described in Appendix A. Assuming that the normalized 



-7- 



variables x and y are of the order of unity and using the scaling relation (|lq), we obtain the final 
results of the expansion up to 0(a^ n ) flChodorowski fc Lokas 1997 , Bernardeau & Kofman 1995): 



A 



- 2 /2 f(5 2 \\ 1/2 



2vr \{Sl 



A,, //,(.<•) h- ^# 2 (s) + ^# 4 (x) 



A11A40 A30A21 



24 



12 



^(*) + ^ 



and 



V(5 n 



-x 2 /2 



2vrv^ 



A 



1 + '-^H 3 (x) + <J '-^-H^x) + ^# 6 (z 



A40 



A 



30 



6 



24 



72 



where H n {x) is the Hermite polynomial defined by 



H n (x) 



x 2 /2 ( °_ 

dx 



-x 2 /2 



(21) 



(22) 



(23) 



Substituting the results (21) and (|^) into equation (|i7|), the conditional mean {5 g )\s m in the 
weakly non-linear regime becomes 



(^)k = V^) [MlX + ^(x)+^ 2 (x)]. 
where T n ix) denotes the terms of the order 0{a'^ n ) given by 

F x {x) = ^H 2 (x) + ^^[H 4 (x)-xH 3 (x)}, 
2 o 

T 2 (x) = ^H 3 (x) + ^^[H 5 (x)-H 2 (x)H 3 (x)] + ^^[H 5 (x)-xH 4 (x)} 



(24) 



72 



H 7 {x) - xH 6 (x) - 2H 3 (x)H i (x) + 2xH 3 (x) 



We should keep in mind that Ajj still contains the terms higher than Ofo^ 3 ) an d we do not 
approximate the variance of galaxy distribution (5 2 ) in deriving the expression (p4|). 

We first consider the conditional mean to the order 0(a m ). Since we can drop the higher 
order correction to Ajj with i + j > 3, we simply replace Xij with S^a^^ 2 . As for the variance 
and the covariance of S m and 5 g , the definition ( pX| ) gives 



(5 



cr r , 



(S, 



- b a m , 



A 



11 



This is correct within the tree-level analysis. Then the weakly non- linear relation (24) up to 
jYi) is written in terms of the variable S m a»s 

(S 9 )\s m = c 1 S m + ^(s 2 m -(S 2 m ))+0(a 2 m ). 

The coefficients c\ and C2 becomes 

ci = br, c 2 = b (S 2 i - rS 30 ). 



(25) 



(26) 
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The expressions (^5|) and ( f26[ ) are our main results. 

To understand the stochastic property in the biasing relation, the mass density-galaxy relation 
should be noted: 



d5 m 8 m P(8 m ,\Sg), 
B 



ns 9 



; B = J d5 m 5 m V(5 m ,5 g ). 



Derivation of non-linear mass density-galaxy relation is similar to that of (<5g)|<j m - For the 
expression (p3|), we only have to replace the variables {S g ), Ay and x with (5^), Ajj and y, 
respectively. We here write down the result in terms of 5 g up to 0(a m ): 



(S m )\ Sg =d 1 5 g + ^- (s 2 - (5 2 g )) + <D(a 2 J 



Then the coefficients d\ and c?2 become 

a r 



b 2 



(S 



12 



rSt 



03 J 



(27) 



(28) 



We observe that the conditional mean (|27|) does not coincide with the result (^) due to the 
presence of the cross correlation r and the fact that Sij is not symmetric. It will be shown in the 
next section that both relations give the same relation when the cross correlation r becomes unity. 

We next clarify the effect of the higher order corrections on the galaxy-mass relation {5 g )\s m - 
When we take into account the ^(x) term, the 0{a 2 n ) contribution requires the third order 
perturbations. In this case, we cannot ignore the O(cr^) correction for Ajj with i + j = 2. Then 
the evaluation of the first term in the right hand side of ( pl|) is modified. Using the definition (|Io|), 
we can write 



(5 2 ) AnX 



br + 0(a 2 m ) 



(29) 



where the order 0(a 2 n ) terms come from the loop corrections: 



< W>> + < W) + (€ 3 4 3) ) - br ( ([4 2) ] 2 ) + 2(*fi)5W) 



Hence we redefine equation ( |29| ) as b ren r ren 5 m to the 0(a 2 n ) terms. The parameters b ren and r ren 
are the renormalized quantities for the biasing parameter b and the cross correlation r given by 
the tree-level perturbation. For Ajj with i + j > 3, the correction terms are higher than 0(af n ), 
which can be verified from the scaling law (|i"6|). Rewriting the expression (|24| ) in terms of S m by 
using the definition (|2"3|), we get the non-linear galaxy- mass density relation in accuracy of 0(a 2 n ): 



(30) 



where the coefficients 5j are 

C\ — brpn.V' 



ren' ren 



S30 5*21 



rS 2 



30 



31 



10 J 
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C2 = -(S^i - rS 3 o), 

Q = 1(5*31 - rS i0 ) - -(£30 £21 - rSf ). 
o 2 

We immediately see the difference between the above result (^) and the tree-level result (|25|). 
In addition to the contribution of 5^, the coefficient c\ has the extra terms except for the 
renormalized parameters, which are absent from the coefficient c\. The result says that the linear 
biasing relation in the weakly non-linear regime is different from the linear theory prediction even 
for S m <C 1. Our results are in good agreement with the analyses of the two point correlation 
function for the galaxy distribution (Scherrer & Weinberg 1998). 



3. Time evolution of galaxy-mass density relation 

The previous section has been devoted to the discussion of the galaxy-mass density relation. 
In this section, we investigate the time evolution of it. The evolution equations and the initial 
conditions are given in section 3A. Within this model, we calculate the time-dependence of the 
joint moments Sij explicitly and the relationship between the conditional means {8 g )\$ m and 
(S m )\s is clarified in Sec |3.2| , In section |3.3| , we will study the time evolution of the conditional 
mean (S g )\s m in order to reveal the influence of initial stochasticity. 



3.1. Evolution equations and initial conditions 



The time evolution of the galaxy biasing induced by the gravity has been studied by several 



authors ( [Fry 199q , [Tegmark fc Peebles 1998| , [Taruya, Koyama fc Soda 1999|) . In these articles, the 
density fluctuation 5 m is assumed to be evolved following the equation of continuity and the Euler 
equation. On large scales, the total mass distribution moves along the irrotational velocity flow. 
It is convenient to define the velocity divergence 9 = V • v/(aH), where a is the scale factor of 
the universe and H is the Hubble parameter. For the galaxy distribution, we consider the epoch 
after galaxy formation and assume that the merging process is negligible. The distribution 5 g 
should satisfy the equation of the continuity whose velocity field is determined by the gravitational 
potential. Then the evolution equations for the total mass and the galaxies become 



05 



1 



+ H + _ V ■ (5 m v) 
at a 



0. 



dd 



n a 

1 ~ 2 + 3H 2 J 



] H + ?-Htt5 m + -4rj 
1 a z H 



V • (v ■ V)« = 0, 



(31) 
(32) 



and 



§ + fffl+-v-(V) = o. 

at a 



(33) 
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The variable A is the cosmological constant and f2 is the density parameter defined by 

n = —4. (34) 

Equations (|3T|) , (|32]) and (|3^) are our basic equations for the time evolution of 8 m and 5 g . The 
same situation has been studied by Fry (1996) in the case of the deterministic biasing (g). 

Next, we consider the initial condition given at an initial time t{. The total mass fluctuation 
5 m is produced during the very early stage of the universe, whose initial distribution may have the 
random Gaussian statistics. We regard such fluctuation as A m (x). The gravitational instability 
induces the deviation from the Gaussian statistics and the galaxy formation does not affect the 
evolution of 5 m on large scales. We give the initial condition 5 m = /(A m ) from the perturbative 
solutions by dropping the decaying mode, which leads to the non-local form of the function y(A m ) 
(see Appendix B). 

On the other hand, the fluctuation of the galaxy number density is induced by the galaxy 
formation. To give the initial non-Gaussian distribution 5 g = g(A g ), we need to know the galaxy 
formation processes. Currently, it is not formidable. Here, we treat g(A g ) as a parameterized 
function, whose unknown parameters are determined by the observation of galaxies. Assuming 
g(A g ) as a local function of A g , we take 

g(A g ) = A g + ^(A 2 g -(A 2 g )) + .... (35) 

Due to the assumptions stated above, in the perturbative regime, the relation between the 
galaxy and the total mass distribution can be characterized completely if the stochastic property 
of the Gaussian variables A m and A g is given. It is expressed in terms of the three parameters : 

-o 2 = <A^>, bl = ^- r = (AgAm) . /2 , (36) 



which is equivalent to giving the matrix G in the generating function (y). The variables bo, tq and 
o"o corresponds to the initial biasing parameter, the initial cross correlation and the initial variance 
of the total mass, respectively. The variance o"o is related to the initial power spectrum P(k) as 

/d?k 2 
-^—^Wl{kR)P(k) ■ W(kR) = jj-^[sm(kR) - kRcos(kR)}, (37) 

where W(kR) is the top-hat window function in the Fourier space. Hereafter, we simply assume 
that the parameters 6q and rg are constant, which is consistent with the fact that there is no 



evidence of the scale-dependent biasing on large scales flMann, Peacock &: Heavens 199 
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3.2. Variance, covariance and joint moments 



Within the above prescription, we will investigate the time evolution of the galaxy-mass 
density relation in the weakly non-linear regime. In the lowest order of the perturbation, the 



conditional means ( P5|) and (27) are written in terms of the variance of the total mass a m , the 
biasing parameter b, the cross correlation r and the joint moments Sij with i + j = 3. 

The parameters a m , b and r are obtained from the linear order solutions of evolution equations 



(31), (32) and (p3|). The solutions satisfying the initial conditions in the previous subsection 
become 



A m (x) D(t), 

A m (x) (D(t)-1)+A g (x), 



where the function D{t) denotes the solution of growing mode by setting D(ti 
(Peebles 1980|) 



(38) 
(39) 

1, which satisfies 



D + 2HD - -H 2 nD = 0. 



(40) 



We have D(t) = a(t)/a(ti) in Einstein-de Sitter universe (O = 1). Substituting (|38|) and (gS|) into 
(|i~0|), the time dependent parameters are evaluated from the knowledge (||) and (36) as follows 
(Fry 1996| , [Taruya, Koyama fc Soda 1999| , |Tegmark fc Peebles 19981) : 



m 

r{t) 



(T D, 



(D_l)2 + 2 6 r (Z)-l) + 62 



b-\t) 



D 

D-1 + b r 
D 



(41) 



We see that the initial cross correlation r$ = 1 leads to r = 1. Thus, in our prescription of the 
time evolution, the stochasticity in the galaxy biasing comes from the initial conditions. 

For the joint moments Sij, we must calculate the second order perturbations by solving the 
evolution equations. In Appendix B, the solutions of second order perturbation are summarized 
and the computation of Sij is explained. The time evolution of the joint moments are parameterized 
by the initial conditions, bo, ro, h, and the power spectrum P(k). The results in Appendix B are 



S30 
S21 




(42) 



D 2 -l 2(D - 1 + bpfp) 1 
7 / I D 2 D J 



(43) 
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S 



12 



3b 2 (t) 
+ (6 - 7) 



34 

y 



7 



(Zj_l + fr orQ )2 | 2(D2_ 1 )( jD _ 1 + fc oro ) | 



D 2 



D 3 



2( J D-l) 2 (6 r -l) , ( J D-l)((6 r -l) 2 + 6 2 -l) 



D 3 



+ 



L> 3 



2b r (bor (D-l)+b 2 ) 



S, 



03 



34 

y ~ 7 



L)3 



(D 2 -l)(D-l + bor ) s 
D 4 



(44) 



+ (6 - 7) 



(D - 1){D - 1 + 6 r ) {(6 r -l)D + b 2 + l- 2b r } 



D 4 



2\ 2 



+ h I w 



{b r (D-l) +b 2 } 
D 4 



(45) 



where the numerical value Iw, which comes from the non-Gaussian initial distribution of the 
galaxies, is defined by 



Iw = o"o 



' 2 Wfl(|fc! + k 2 \R)W R {k x R)W{k 2 R)P[k x )P{k 2 ) 



It has been checked by the Monte Carlo integration that Iw is almost equal to unity. Hence, we 
can regard h as the initial skewness of the galaxy distribution. The variable 7 is given by 

d 



7 



[loga 2 (R)]. 



d(logR) 

For the power spectrum with the single power-law P{k) ocP, we have 7 = n + 3. 

By definitions (|l0|) and (|l6|), the moment 6*30 is identical with the skewness of the total mass 



5*3,™ (see |15[). The skewness of the galaxy distribution is related to the moment S03 as 
Ss t9 = Sos/b(t). In the case of tq = 1, using these facts, and equations fl4^)-(|45|), we can get the 
relation 

c 2 = 6(521 - rS 30 ) = -(65 3iS - 5 3 , m ). 

Thus, in the deterministic case, the coefficients c\ and c 2 can coincide with 6 and 62 given by 
equations (J2J) and fli~3|). It should be recognized that the galaxy-mass density relation defined by 
equation (0) does not coincide with the conditional mean (5 g )\s m in general. 



Using the results (42)- (45), we have another important conclusion for the galaxy-mass density 
relation. We can obtain 



S21 — S; 



30 



12 



S, 



03 J 



for ro = 1. Taking 5 9 = 6<5 m , this means that the approximation of <5 m by the perturbative 
inversion of {8 g )\s m coincides with the conditional mean {6 m )\$ (eq.^7|). However, if we take into 
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account the higher order corrections, we can expect that the conditional means {S g )\s m and {5 m )\$ g 
does not become equivalent even in tq = 1 case, because of the higher order contribution to the 
linear order coefficient of the conditional means (see eq.|3(| below). 



3.3. Evolution of (S g )\ Sm 

Now consider the time evolution of the conditional mean (S g )\s m to see the influence of the 
initial stochasticity on the galaxy-mass density relation. For brevity, we only describe the analysis 
in the Einstein-de Sitter universe (fio = 1). 

Using the results (PH), (42) and ([43]), the time evolution of the coefficients c\ and c<i can be 
examined. In Fig.l, we plot the coefficients c\ and ci given by equation (^) as a function of 
the expansion factor a(t). As the initial parameters, we chose bo = 2.0, tq = 0.8 and h = 3.0 at 
a(ti) = 1 in Fig. la. The solid line shows the coefficient c\. The long-dashed, the short-dashed and 
the dotted lines denote the time evolution of the coefficient C2 with the spectral index n = —2, 
—1.5 and —1 for the single power-law P(k) oc k n , respectively. The parameters in Fig. lb is the 
same as Fig. la, except for the initial cross correlation for which we take tq = 0.1. These figures 
show that the coefficients c\ and C2 approach to unity and zero respectively, independent of 
the initial parameters and the spectral index. This asymptotic behavior can be ascribed to the 
attractivity of the gravitational force (Fry 1996, Tcgmark & Peebles 1998, [Taruya, Koyama fc Soda 



1999). Fig.lc shows the illustrative example with the initial parameters &o = 4.63, tq = 0.2 and 
h = 6.96 at a{ti) = 1. We set the spectral index n = —1.41. If we identify the initial time a = 1 
with the redshift parameter z = 3 assuming the Einstein-de Sitter universe, the skewness and 
the bi-spectrum at present time (which corresponds to a = 4 in our case) provide the consistent 



results with the observation of the Lick catalog ( Taruya, Koyama fc Soda 1999 ). Although the 



initial value of the coefficient C2 is comparable to that of c\ because of the large initial skewness 
h, it rapidly decreases and becomes negligible due to the small initial cross correlation r$. As a 
demonstration, we evaluate the galaxy-mass density relation (8 g )\s m in the (5 g , <5 m )-plane with the 
same initial parameters as Fig.lc. The galaxy-mass relation is plotted in Fig. 2 by choosing the 
initial variance of the total mass do = 0.1. Each line in Fig. 2 represents the snapshot at a = 1 
(solid line), a = 2 (long-dashed line), a = 4 (short-dashed line) and a = 8 (dotted line), respectively. 

Fig.l and Fig. 2 say that the coefficient C2 is usually smaller than the linear coefficient c%. We 
can also confirm this fact for various initial parameters. This result means that the non-linearity 
in the conditional mean is negligible on large scales. Therefore the result indicates that the 
relation between the galaxy and the total mass distribution is well-approximated by the linear 
relation 5 g = brb m . If we apply this fact to the estimation of the higher order moments for galaxy 
distribution, we expect that the moments of the galaxy are simply related to those of the total 
mass multiplied by the factor inferred from the linear relation. That is, regarding (S g )\s m as S g , 
we obtain 

<(*,)»> ~ (brT((5 m T). (46) 
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However, we should keep in mind that there exists the scatter in the galaxy-mass relation. In 
section 2.2, we have defined the biasing scatter e (see eq.|12|). In the tree-level analysis, the 



variance of the biasing scatter becomes 



<1> =(1-^ + 0(0. (47) 



The variance vanishes only when tq = 1. We can also calculate the third order moment. The 
result up to the second order perturbation is 



< 



SQ3-r- i S 30 -3r(S 1 2-rS 2 i) +0(a 2 m ), (48) 



where we used c\ = br. Substituting the expressions (fi^)-(H^) and ([II]) into equation (jig), it is 
easily shown that the third order moment of the biasing scatter also vanishes in the deterministic 
case, ro = 1. From equations (f47|) and ([48]), the non-Gaussian scatter is expected to affect the 
simple relation of the higher order moments ([46]) . 



To see the influence of the biasing scatter on the evaluation of higher order moments, we 
examine the skewness of the galaxy distribution. In Fig. 3, we plot the time evolution of S% ig as a 
function of the expansion factor. For each figure, the solid line is the correct skewness obtained 
from the joint moment divided by b{t) and the dashed line represents the skewness S^ t g 
deduced from the galaxy-mass density relation {5 g )\$ m . §3 tg can be evaluated by equating the 
conditional mean {5 g )\$ m with 5 g . The non-linear relation between 5 m and 5 g ( |25| ) implies 

S3, 9 = -(S3,m + S-), (49) 
Cl Cl 



which is valid within the tree- level approximation ( Fry Gaztahaga 1993] ). We have already 
obtained the skewness of the total mass S^m = %- — 7 and found that 5"3 im is usually larger than 
the coefficient 3(c2/ci). Fig. 3a and Fig. 3b have the same initial parameters as Fig. la and Fig. lb, 
respectively, except for the power spectrum specified as the index n = —2. As we know from the 
previous subsection, S^g exactly coincides with the skewness 53 i9 only in the deterministic case, 
i.e, ro = 1. The figures show that 6*3^ differs from £3^ in the presence of stochasticity clearly. For 
the smaller initial cross correlation ro, the deviation of £3^ from Ss tg becomes more significant 
(Fig. 3b). Because of the small contribution of the non-linear coefficient C2, £3^ approaches to 
the skewness of the total mass S3 : m = 3.86 more rapidly than S3 ;9 (see eq. f49|). For the same 
parameterization as used in Fig.lc, Fig. 3c shows that £3^ at present (a = 4) underestimates the 
correct value of the skewness for the galaxy distribution due to the rapid relaxation to S3 j7n . At 
o = 4, the skewness 6*3^ from Fig. 3c becomes 4.55, while we have S3 t9 = 3.50. The difference 
cannot be neglected. Thus there exist the cases that the galaxy-mass density relation as a 
conditional mean leads to incorrect result. This feature is common to the stochastic biasing. 
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4. Conclusion 



In this paper, we have analytically studied the galaxy-mass density relation in the framework 
of the stochastic biasing. Under the assumptions (|5|) and (]6|), in the weakly non-linear regime, we 
derived a general formula for the galaxy-mass density relation as a conditional mean. The main 



result in our analysis is equations (25) and (|2q). We have seen that the higher order contribution 
to the weakly non-linear galaxy-mass density relation can shift the coefficient of the linear 
galaxy-mass density relation, in addition to the non-linear term proportional to 5^ (see eq. |30fl ) . 
This agrees with the analyses of the two point correlation function ( [Scherrer fc Weinberg 1998[ ). 
Using the formula for the conditional mean, we have further investigated the time evolution of 
the galaxy-mass density relation. To develop the analysis, we have made the assumptions: (i) the 
galaxy formation and the merging process can be ignored; (ii) the initial distribution of galaxies is 
given by the local function (eq.|35]]), while the total mass has the non-local initial condition. Our 
conclusions can be summarized as follows: 



• The conditional mean could be different from the deterministic biasing relation (compare the 
relation [Q and P6[ ). In the presence of the scatter, the perturbative inversion of (S g }\s m 
does not recover the conditional mean {6 m )\s ■ 

• The time evolution of the conditional mean (5 g )\$ m shows that the non-linear term of the 
conditional means usually becomes negligible. This suggests that the galaxy-mass density 
relation could be approximated by the linear relation. 

• We have found that the time evolution of the skewness deduced from the conditional mean 
(6 g )\s m differs from that of the correct skewness S^^g in the case of the small cross correlation. 
This indicates that the stochasticity could have an important role in the estimation of the 
higher order moment of galaxies using the conditional mean. 



The assumptions (i) and (ii) might not be valid in more realistic situation. The merging 
process might become important in a real universe and the galaxy formation would introduce the 
non-local initial condition. Nevertheless, we believe that our qualitative conclusions will not be 
altered even if these processes are taken into account. It would be interesting to incorporate these 
effects into the time evolution of the stochastic biasing. We will attack this issue in the future. 

The conditional mean (S g )\s m is an important quantity to construct the analytic biasing 
model which is determined by the clustering properties of the halos ( [Mo fc White 1996| , Mo, Jing 
fc White 1997, patelan, Matarrese fc Porciani 1998 , 3heth & Lemson 1998). In the presence 
of the stochasticity, we must treat the conditional mean carefully when comparing it with the 
observation of the galaxy statistics, and vice versa. That is, the linear and non-linear biasing 
parameter (6, 62) estimated from the observation of the skewness or kurtosis does not necessarily 
give the mean biasing relation. They might provide a signal for the stochastic biasing. In the 
previous paper ( Taruya, Koyama fc Soda 1999[) , we found that the parameter 62 deduced from 
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the skewness can become negative in the presence of the scatter. The recent observation from the 
Southern Sky Redshift Survey shows that the non-linear biasing parameter 62 estimated by the 
skewness may become negative for the biasing parameter b > 1 ( |Benoist et al. 1998| ). This might 
give an observational evidence for the stochastic biasing as long as our prescription is correct. 

The stochasticity is problematic when we get the relation between the galaxies and the total 
mass from the observation. As discussed by Dekel & Lahav (1998), the situation may become 
more complicated in the redshift space. This makes it difficult to determine the cosmological 
parameter from the observation of the velocity field such as POTENT ( |Bertschinger Dekel 



1989|) . For the analysis in the redshift space, Scoccimarro, Couchman & Frieman (1998) studied 
the influence of the non-linear biasing on the bi-spectrum of the galaxies in the deterministic case. 
We must investigate how the stochasticity and the non-linearity affect the galaxy biasing in the 
redshift space. In particular, the higher order statistics such as the skewness and the bi-spectrum 
should be explored. Extension of our formalism to the redshift space is straightforward and the 
analysis is now going on flTaruya fe Soda 1999| ). 

Alternative approach to understand the stochastic property of the galaxy biasing is to seek 
the physical origin of the stochastic biasing itself. Blanton et al. (1998) explored the hidden 
variable to reduce the stochasticity in the relation between galaxies and total mass using the 
hydrodynamical simulation. They found that the scatter around the conditional mean of the 
galaxies becomes small if the temperature dependence is taken into account. To combine this 
approach with ours is also interesting. 



This work is partially supported by Monbusho Grant-in- Aid for Scientific Research 10740118. 
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Appendix A: Calculation of Edgeworth series 

In this appendix, we explain the calculation for the derivation of the galaxy-mass density 
relation conditional mean in section 12.31. 

We start to write down the generating function W expanded in powers of a m . From the 
scaling relation (([(]), the generating function up to the order 0(a 2 n ) is obtained in terms of the 
normalized variables (|i~8|) : 

W(/i, u) = -^(^ 2 + 2Au/ii/ + v 2 ) 



+ 



+ 



g {AsoM 3 + 3A 2 i(i/i) 2 (ii/) + 3X 12 {ifM){iu) 2 + A 03 (^) 3 } (1) 
jt {a 4 oM 4 + 4A 3 iM 3 (i^) + 6A 22 (^) 2 (^) 2 + 4X l3 (ifx)(iu) 3 + A 04 (^) 4 } 



24 

+ 0{al). 

Note that the first line in the right hand side of ([[]) is of the order of unity. The second and 
the third lines are the terms of 0(a m ) and 0(a 2 n ), respectively. Using the expression (|l|), the 
conditional mean {S g )\s m given by equation ( |l7| ) is calculated in the following way. For the 
numerator A, we rewrites the expression (|l^) to 



A 



ov 



Integrating by part and substituting the expansion (|]) into equation @, we get 

1/2 r h 



A 



27T I (51 



I 

2^ 



1 + 



+ 0(a. 



Keeping the terms up to 0(a 2 n ), this becomes 



A=± 



1/2 



27T I <(& 



An(-^j+^(~ 



+ 



A11A40 A30A21 



24 



12 



d;i; 



^ \ , ^11^30 



dx 



+ 



An A 



da: / 



+ 



72 



da: / 



11^30 

6 

+ 0{al) 



dx 



A31 
6 



_d^ ! 

dx, 



dfj,e 



where we replaced (i/x) n with (^—-^j • Performing the integration of \i and the differentiation 
with respect to x, we obtain the expression (21). 
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The calculation of the denominator in equation ([[7]) is similar to that of the numerator. 
Substitution of the expansion ([[]) into equation ( p0| ) yields 



V{5 n 



d\x e 



-i^,x—y?/2 



1 i ^30/. \3 , 

i + -g-M + 



Repeating the same manipulation as in the above, we have 



V(5„ 



1 



6 



d 



dx 



+ 



A40 
24 



d 



dx 



+ 



A 2 



30 



72 



dx, 



-x 2 /2 



Using the definition of the Hermite polynomials (|23|), we obtain the expression (p^) 



Appendix B: Second order perturbation and joint moments 



When we investigate the time evolution of the conditional mean {5 g )\s m and (5 m )\s g in the 
tree-level analysis, in addition to the biasing parameter b and the cross correlation r, we need the 
joint moments Sij with i + j = 3: 



S: 



30 



(I)l2x(2)\ 
m I Cm , 



([5m 



(1) 



s- 



2\2 



21 



([5 



«] 2 5f)+2(^^4 1) ) 



6(t) <[4 



5 



12 



(4 2) [4 1) ] 2 ) + 2(^ ) 4 1) 4 2) ) 



5, 



& 2 (i) ([^ } ] 2 } 2 ' ~° 3 ([S^] 2 ) 2 ' 
If we evaluate these quantities for the smoothed density fields, the integration including the 
window function Wr(x) must be computed. Although the integration in real space is difficult for 
the higher order moments, it is known that the calculation is tractable in the Fourier space because 
of the useful formulae ( |Bernardeau 1994] ). Using these formulae, we shall write down the solutions 
of second order perturbation in the Fourier space. The variables 5 m ,S g and 6 are expanded as 



(1) 



d k ~ . . ikx /1 / \ 

(27I") 3 > 9 ' ' ' e ' v 35 ' ) 



d 3 k s 



(2tt)5 



9(k,t)e 



ikx 



Then the second order solutions satisfying the initial conditions are obtained (Fry 1984, Taruya, 
Koyama fe Soda 1999j ): 



d 3 k' 



D 2 (t) %R,(k', k-k') + -K{k - k', k') - -C(k', k - k' 



+-E{t)C(k',k-k') 



A m (k-k')A m (k'), 



(2) 



8^(k,t) 



d 3 k' 



r d°k 

+ (D(t) -I) J 7^3 W, k - k')(A m (k')A g (k - k') - A m {k')A m {k - k')) 



+ 



h 
6 



d 3 k' 
(2vr) 3 



(27T) 

A g (k')A g (k-k')-(A 2 g )5 D (k') 



(3) 
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where 

v( U u\ i , ( fc l- fc 2) r(u u\ 1 ( fcl -fc2 ) 2 
|fci| |fei| z |fc2r 

The solutions (|2|) and (|3]) contain the function E(t) which satisfies E{ti) = 1. This is the 
inhomogeneous solution of the following equation: 

E + 2HE - -H 2 nE = 3H 2 QD 2 + -D 2 . 

2 3 

In Einstein-de Sitter universe, we have 

E(t) = ^D 2 (t). 

It is known that the O and A dependences of E/D 2 is extremely weak ( [Bernardeau 1994] ). 
Therefore, we proceed to evaluate Sij by replacing E(t) with (34/21)D 2 (t). Substituting the 
solutions @ and (|3|) into the expression (|l|) and using the formulae below, we get the results 



K(k u k- 2 )P{ki)P(k 2 )W R (\k l + k 2 \R)W(k 1 R)W{k 2 R) = o%{R) ( 1 



where 



(2vr) 3 (2vr) 3 ' 

^^£(fci,fc 2 )P(fci)P(fc2)Wji(|fci + fc2|lJ)W(fci22)W(fe/2) = 3^(P), 



^ - 2 



d(logP) 



[log^(P)]. 



Wn{kR) and o"o are the window function in the Fourier space and the initial variance of the total 
mass defined by equation (|37|), respectively. 
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Figure Caption 

Fig.l The time evolution of the coefficients c\ and c% is evaluated as a function of the expansion 
factor a(t) in the case of the Einstein-de Sitter universe. For each figure, the solid line 
shows the coefficient c\. The long-dashed, the short-dashed and the dotted lines in Fig. la 
and lb denote the time evolution of C2 with the spectral index n = —2, —1.5 and —1, 
respectively. The dashed line in Fig.lc also shows the coefficient ci, but we set the spectral 
index n = —1.41. The initial parameters for each figure are as follows: (a) bo = 2.0, r$ = 0.8 
and h = 3.0; (b) b = 2.0, r = 0.1 and h = 3.0; (c) b = 4.63, r = 0.2 and h = 6.96. Note 
that the parameters in Fig.lc gives the same skewness and the bi-spectrum as obtained from 
the Lick catalog data at a = 4 if the universe is assumed to be the Einstein-de Sitter universe 
(see section j3lj|) . 

Fig. 2 The time evolution of the galaxy-mass density relation {6 g )\$ m is plotted in (S g , 6 m )-space. 
The initial parameters are the same as Fig.lc. Each line represents the snapshot evaluated 
at a = 1 (solid line), a = 2 (long-dashed line), a = 4 (short-dashed line) and a = 8 (dotted 
line). 

Fig. 3 The time evolution of the skewness for the galaxy distribution is shown as a function of the 
expansion factor a in the Einstein-de Sitter universe. The solid line is the joint moment S03 
multiplied by the factor l/b(t), which is identical to the correct value of the skewness S^^g. 
The dashed line denotes the skewness £3^ deduced from the galaxy-mass density relation 
(|25|), which is evaluated from equation (|49"1), Fig. 3a and 3b respectively have the same initial 
parameters as Fig. la and Fig. lb, except for the spectral index specified as n = —2. As 
for Fig. 3c, the initial parameters are the same parameters as Fig.lc, which gives the same 
results as obtained from the Lick catalog data if we identify the initial time a(t) = 1 with 
the redshift parameter z = 3 in the Einstein-de Sitter universe. 




Fig. la 




Fig. lb 
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